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Abstract. The Rayleigh-Modal method is used to calculate the electromagnetic field 
within the grooves of a perfectly conducting, rectangular-shaped ID diffraction grating. 
An enhancement coefficient (77) is introduced in order to quantify such an energy 
concentration. Accordingly, rj >1 means that the amount of electromagnetic energy 
present within the grooves is larger than that one will have, over the same volume, if 
the diffraction grating is replaced by a perfectly reflecting mirror. The results in this 
paper show that 77 can be as large as several decades at certain, often narrow, ranges of 
wavelengths. However, it reduces to approximately 20% under sunlight illumination. In 
this latter case, such values are achieved when the optical spacing between the grooves 
dn is greater than 500 nm, where d is the groove spacing and n is the refractive index 
of the substance within the grooves. For dn smaller than 500 nm the enhancement 
coefficient turns negligibly small. 



PACS numbers: 2.79.Dj, 88.40.H-, 88.40.fc, 42.79.Pw, 42.79. Qx 
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1. Introduction 

It is well documented that electromagnetic fields within the grooves of diffraction grating 
can be largely increased when illuminated by light [H [21 El H]. This phenomenon 
was recently proposed as way to enhance the absorption of light in PV cells and 
optoelectronic devices [HI El [71 [8] . In fact, for the purpose of increasing the efficiency 
of photo-cells where sunlight may not be absorbed so easily, the field enhancement 
appears to be a useful approach. In this realm, however, rather than differential 
one needs integrated figures. This is so because large field enhancements, as those 
previously reported, are not enough if not accompanied by a net increase of the light 
energy integrated over a representative volume of the cell, along the pertinent range of 
wavelengths, incidence angles and states of polarization. 

In a previous paper of the Authors [9], the field enhancement was analysed for 
the more general case of beams of monochromatic light arriving to the grating along 
several directions and different states of polarization. There, an enhancement coefficient, 
namely 77, is introduced, so that rj >1 means that the amount of electromagnetic energy 
present within the grooves is larger than that one will have, over the same volume, if the 
diffraction grating is replaced by a perfectly reflecting mirror. In this paper, however, 
the previous study is extended to the case of illuminating the grating with solar light. 
The results of these calculations show that, when using monochromatic, polarized and 
well-collimated beams of light, rj exhibits a series of peaks at certain wavelengths, where 
the enhancement coefficient can be as large as several decades, or even grater. However, 
after taking an average over incidence angles and polarization states, rj is significantly 
reduced as it may reach then values which can hardly be larger than approximately 
3. If the previous results are also averaged over the solar spectrum, the enhancement 
coefficient is further reduced to approximately 1.20. 

Although these results are to some extent discouraging, it does not mean that 
the diffraction grating structure cannot be advantageously used in designing photo- 
electronic devices. This is so because, depending on the case, the aforementioned 20% 
gain may suffice. In this regard, the present study could be interesting to those who may 
need increase the absorption of light within a solar cell and might possibly be planning to 
incorporate a diffraction grating for such a purpose. This paper is organized as follows: 
the basics equations necessary to obtain the field within the grooves of a diffraction 
grating are derived in Sectional The results of numerically calculating the enhancement 
coefficient as a function of the various parameters in the model are shown and discussed 
in Section [3] and, finally, Section [4] contains a summary and concluding remarks. 

2. Basic equations 

The calculation model is sketched in figure [U There, one can observe the incident light, 
a linearly polarized plane-wave, arriving to the grating from above at an angle 9i with 
respect to the ?/-axis and azimuthal angle (fi with respect to the x-axis. The light has a 
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Figure 1. A linearly polarized electromagnetic wave, with wave vector kj, electric 
field Ej, incidence angle 9i and azimuthal angle (fi, arrives from the upper half-space 
upon a rectangular-shaped diffraction grating. Grooves have a spatial period d, width 
b and depth h, and cj) denotes the polarization angle with respect to the (x-z)-plane. 



wavevector kj = (ki X , ki y , ki Z ), which is related to the wavelength in the upper medium 
(A) as |kj| = ki = 2n/\. The electric field amplitude Ej is assumed to be perpendicular 
to kj and <fi is the angle formed by Ej and the (y = 0)-plane. 

The diffraction grating is assumed to be an infinite set of equally spaced, 
rectangular-shaped grooves made on the surface of a perfect conductor. The spatial 
period of the grooves is d, whereas the width and depth are b and h, respectively. 
Finally, the upper medium is assumed to be a non-lossy dielectric material that fills the 
upper part of the grating, up to a front surface which will not be taken into account in 
this calculations. 

The way the electromagnetic fields within the grooves are calculated was already 
shown in Ref. [TJ |9], however, for the purpose of the present analysis, it is important to 
re-do it once more. 

In the first place, the Rayleigh expansion is assumed for upper space, 

E (y>0) ( r ) = E.g^-r + R n e iq - r for y > , (1) 

n 

where q n is the wave vector of the ra-th order reflected beam and R n = (R nx , Rnyi Rnz) 
is the corresponding electric field vector. Similarly, q n = (q n x,Qny,k iz ), where q xn = 
k ix + 27rn/d, n = 0, ±1, ±2, . . ., and q ny = +\Jkf — q\ x — kf z . Notice that the plus sign 
in front of the square root means that the positive, either real or imaginary solution to 
the square root has to be used. 



Light trapping within diffraction gratings 

Within the grooves one has the so-called modal waves, namely, 



4 



^ <0) (r) = '£ £ ^me i ^ z+Nk ^X' m>N (x)Y m (y) 

m,N 

E i V<0) ( r ) = E £ v>™ e^ z+Nk ^ X m , N (x) Yi(y) (2) 

m,N 

Ely <0 \r)=tJ2^rne tik ~ z+Nk - x) X m , N (x)YUy) 



m,N 



where 



X m N(x) 



sin [K mx (x — X]\f)] for xn < x < x n + b 
otherwise 



/ J cos [K mx (x - x N )} for x N < x < x N + b 

X ^ = { o otherwise (3) 

Y m (y) = sin [K my (y + h)\ and Y,' m (y) = cos [K my (y + h)\ , 

where K mx = mir/b for m = 0, 1, 2, . . . and K my = + ^kf — kf z — K^ x , where, again, 
the positive solution to the square root must be taken. Similarly, a; at denotes the in- 
coordinate of the right wall in the N-th groove, i.e. x^ = Nd for iV = 0, ±1, ±2, . . .. 

Both, the continuity of the x and z component of the electric field along the groove 
aperture, and the boundary conditions of electromagnetic fields on the surface of a 
perfect conductor (see Refs. fTUl [TT] )- allow one to write 

f+oo 
£ x,m cos (K mx x) sin (K my h) 0<x<b 
m=0 

b<x <d 

(4) 



and, 



+oo 

d e iq n *x = \ i Yl £z > m sin ( K ^x) sin (K my h) < x < b 

m=0 



E- e ikixX + R 



n= -°° { b<x <d 

(5) 

By taking the Fourier transform of these equations in the x variable and after rearranging 
terms, one obtains 



Rnx E{ x 5 n Q ~\- — ^ ^ £x,m Cm,n SU1 (K m y}l) (6) 

m=0 

i 

Rnz — —Ei z 5 n fi + — £-z,m £m,n SU1 (K my tl) (7) 



m=0 
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where, 



dx e iqnxX sin (K mx x) and C m , n = dx e tQnxX cos (K mx x) . 
n Jo 



Similarly, the continuity of the x and y components of the magnetic field along the 
grooves aperture yields 



-» \^i y + kf z ) E iz + ki X ki Z E ix 
2 — y - : exp [ik ix x) 



and, 



where 
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It must be noticed that no equations for the y-components of the field are necessary 
since the null divergence condition, i.e. V • E = 0, implies that only two, out of the 
three components of both E^ >0 - ) and E^ <0 ' ) are independent. Finally, the values of £ x>m 
and £ 2 m can be found by evaluating Eqs. fl9 HTU|) over a finite set of equally spaced x's 
along the groove aperture, and the resulting system of linear algebraic equations is then 
solved by resorting to the Gauss elimination algorithm in Ref.[12]. Once £ x>m and £ z , m 
are obtained, the reflection amplitudes, i.e. R nx and R nz , can readily calculated from 
Eqs.OMZD. 
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Having arrived at this point, the enhancement coefficient (77) is introduced, namely, 



7/ 



2A|E,-P 



alx / dy\E(x,y)\' 



(12) 



'0 J-h 

where |Ej| is the electric field amplitude of the incoming light, E(x, y) is the electric field 
within the groove and A is the cross sectional area of the groove, i.e. A = hb. Notice 
that no integration over the z-coordinate is necessary since, for ID gratings |E(x, y)\ 2 
does not depend on z. 

As was previously mentioned [9], such a coefficient allows one not only to quantify 
the field enhancement but also, and most important, it tells when there is a net gain 
of electromagnetic energy within the grooves. Actually, one can readily see that 77 = 
1 denotes the case for which, on an average, the electrical energy density within the 
grooves is the same as that present above, far away from the diffraction grating. 
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Figure 2. Enhancement coefficient for a h/b=l (O)i 2 (A), 4 (y) and 8 (0) diffraction 
grating, illuminated by non-polarized light, uniformly distributed over a 30-degrcc- 
wide, vertical cone. 

It must stressed that the main goal in this paper is the calculation of rj as well as 
analysing its sensitivity upon the various parameters of the model. Before going into the 
results, however, it must be noticed that for the purpose of exploiting the trapping of 
light within the grooves, one can reasonably assume that d = b. By doing so, the number 
of variables is reduced by one and b can be used as the unit of length. Accordingly, apart 
from 9i, (fi, |Ej| and (ft, the solutions to Eqs.fl MTO]) must be functions of bki (or A/6), and 
h/b. The results of numerically calculating these equations under the aforementioned 
assumptions are produced in the following section. It must be mentioned however, that 
all along the following sections, when referring to wavelength it is assumed to be that 
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in vacuum. Only within the numerical code, the wavelength is translated into that of 
Silicon using the refractive index from Ref. [T3] . 

3. Results and discussions 

The enhancement coefficient is calculated for a 30-degree conical beam of non-polarized 
light. The results are plotted in figure [21 where rj appears as a function of the wavelength 
and for grooves of depths h/b—1, 2, 4 and 8. It must be mentioned that these results 
are obtained by calculating the mean- value of r\ in Eq.( fl2l) over the pertinent range of 
polarization states, incidence angles and wavelengths. This is performed by means of a 
Monte-Carlo integration scheme, where relative uncertainties of the order of ten percent, 
or less, is normally achieved. 

As one can readily see in figure EJ the enhancement coefficient exhibits a single 
large peak around X/b ~ 1.8, becoming as large as approximately 3 for the four groove 
depths analysed in this paper. For \/b <1.8, r\ decreases and becomes slightly greater 
than unity for X/b < 1.2. However, one can see that rj falls below unity as soon 
as X/b >2. Curiously enough, the enhancement coefficient does not exhibit a strong 
dependence with the groove depth h and, as a matter of fact, differences between the 
results corresponding to the different groove depths are masked by the noisy aspect of 
the curves. This noise stems in part from the statistical fluctuations in the Monte-Carlo 
calculations, and also from the narrow resonances occurring within the groove which, 
as will be seen below, have been considerably reduced as a result of taking averaged 
values. 

In order to analyse the results in figure [2j rj is calculated for light arriving to 
the grating along the normal direction and for two limiting cases of polarization 
states, namely, those with the electric field perpendicular and parallel to the grooves, 
respectively. The results are plotted in figure [3l where, in order to avoid a busy plot, 
only the results corresponding to h/b=l and 8 appear. 

The results in figure |3]^a) clearly show that when the electric field is perpendicular 
to the grooves, the field enhancements are oscillating functions of X/b, becoming scarcely 
larger than unity over nearly the entire range of wavelengths calculated in the present 
paper. On the contrary, when the electric field is directed along the z-axis, as shown in 
figure EJ^b), 7] exhibits a series of large peaks, regularly spaced, with an amplitude which 
appears to increase with increasing the wavelength. This is so, however, for X/b <2, 
because for X/b >2 the field enhancement falls below unity, and it does at a rate that 
seems to be an increasing function of the groove depth. This however is not at all 
unexpected, since E z must be zero all over the surface of the grating, therefore, a E z - 
polarized light cannot penetrate within the groove as soon as its wavelength becomes 
comparable to, or larger than the groove width. 

Figure H] shows four limiting cases of beams arriving to the grating with 30 degree 
incidence angle These comprise the TE-and TM-polarization state and, and 90 
degree azimuthal angles (ft). The first conclusion one may extract from these curves 
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Figure 3. Enhancement coefficient for h/b=l (full line) and 8 (short-dashed line) 
diffraction gratings, illuminated by a normally incident light with electric field directed 
along the (a) x- and (b) z-axis, respectively. 




wavelength (A/b) 



Figure 4. Enhancement coefficient for h/b=l (full line) and 8 (short-dashed line) 
diffraction gratings. Light arrives at 30 degree incidence angle (0*) upon the grating, 
whereas the azimuthal angle (tpi) and the electric field direction angle (0), both in 
degrees, are: (a) 0, 90; (b) 0, 0; (c) 90, 90; and (d) 90, 0. (see figure Q] for an 
explanation of these angles) 
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is that if the electric field of the light lies on a plane that is parallel to the grooves, 
such as the (b) and (c) cases, then, there appears a cut-off wavelength above which the 
field enhancement falls down below unity quite rapidly. However, if the electric field 
is perpendicular to the grooves, r\ remains greater the unity over the entire range of 
wavelengths analysed in this paper. This is obvious, since no large electric fields can be 
developed within the grooves if they are parallel to the lateral and bottom surfaces. 

As was already observed elsewhere |U [21 E E] , the results in figures H](a-d) show a 
number of peaks, which can reach values as large as several decades. These peaks are 
observed to occur at wavelengths which depend on the groove depth h, the incidence 
and azimuth angles, and the state of polarization. Curiously though, such peaks are 
often so narrow that many of them nearly disappear after taking an average, as seen in 
figure 1. 

Finally, one may find the enhancement coefficient of the grating exposed to sunlight. 
In this case, one must take an average of n over all polarisation states and wavelengths 
in the solar spectrum, 



where A is the wavelength of the light in vacuum, $ S (A ) is the spectral flux density of 
the sun light, r\ is the enhancement coefficient for a given wavelength and after averaging 
over all polarisation states, and Ao m is the largest wavelength of the light which can 
promote electrons from the valence to the conduction band. Notice that the factor Ao 
appearing in both integrals is required in order to obtain the photon density spectrum 
from $e(Ao), which is obtained from the so-called Reference AM 1.5 spectra in [14J. 
It must be also noticed that, since the wavelength in vacuum is used, the dimensions 
of the grating must be scaled using the refractive index n of the medium that fills the 
groove. For the sake of simplicity, such an index is assumed to be constant along the 
solar spectrum and absorption is ignored, as a consequence r\ sun will be a function of 
both the optical spacing of the grooves, i.e. dn, and the aspect ratio h/b. 

Equation (fl3l) is calculated using the Monte Carlo method, and the results are 
plotted in figure |5j There, one can see that the enhancement coefficient becomes larger 
than unity for dn approximately grater than 500 nm, whereas for nd smaller than this 
value, r] sun drops down fairly fast. For optical spacing between the grooves greater than 
500 nm the enhancement coefficient becomes as large as 1.2 and stays around this value 
up to dn=3600 nm, which is the largest <in-value calculated in this paper. Remarkable 
though, in a similar fashion as was previously observed in figure [21 r\ sun appears to nearly 
not to depend on the aspect ratio h/d. 
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Figure 5. Enhancement factor of the diffraction grating under AM f.5 solar 
irrandiance assuming normal incidence. Calculations are performed for several values 
of the optical spacing of the grooves, i.e. dn, and two aspect ratios h/d=l (O) i an d 
8(0).) 



4. Summary and concluding remarks 

The performance of a diffraction grating as a light trapping structure for PV cells 
applications is analysed. To this end, the electric fields produced within the grooves of 
a perfectly conducting, diffraction grating by a well collimated, monochromatic linearly- 
polarized beam of light are calculated. The grating has a (i-period and is assumed to be 
made of infinitely long rectangular-shaped, h-deep and 6-wide grooves. Although, for 
the purpose in this paper it is assumed that d ~ b. Furthermore, in order to properly 
assess the grating, an enhancement coefficient (rj) is introduced, rj is defined as the 
average electromagnetic-energy within the volume of the grooves relative to that one 
will have if grating is replaced by a flat perfect reflector. In this regard, a r] greater 
than unity implies that light is being trapped within the grooves of the diffraction 
grating. Results in this paper agree with previous calculations [TJ |2j |3j HJ U\ that, for 
certain polarization state, wavelength and incidence angle, rj can be substantially larger 
than unity. When using conical, non-polarized beams, such enhancements however are 
reduced and rj can hardly be larger than three. Moreover, these values are observed 
for wavelength in the propagating medium within the range b < A < 2b. For A < b, 
rj appears to be slightly greater than unity irrespective of the polarization state and 
incidence angle, provided this later is not greater than 30 deg. For A > 2b however, the 
enhancement coefficient seems to be always smaller then unity. This is particularly so, 
when the electric field of the incoming light is parallel to the groove direction. Finally, 
the results of calculating the enhancement coefficient of the diffraction grating under 
the so-called Reference AMI. 5 solar spectrum, namely i] sun , show that r] sun can be, in 
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the best of the cases, of the order of 1.2 . This maximum occurs for optical spacing of 
the groove dn approximately equal to 500 nm, whereas for dn < 500 nm rj sun becomes 
negligible small. In all the cases, i] sun does not seem to depend on the aspect ratio h/b. 
It is worth mentioning that, although the structure of the ideal cell sketched in this 
paper is different to that in Ref . [8] , the figures reported by these Authors are similar to 
those in the present calculations. In short, diffraction grating may conceivably act as a 
light trapping structure and, consequently, increase the efficiency of a PV cell, figures 
appear to be around 20% for solar spectrum and groove spacing is expected to be in 
the sub-micron scale. 
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